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The amplitude-probability distribution function of atmospheric radio noise can be 
predicted with reasonable accuracy for a given bandwidth using only the first two moments 
of the noise measured at that bandwidth. This paper presents a method for predicting this 
distribution function for any specified bandwidth from the moments of the noise measured 
at a particuhir bandwidth. 



1. Introduction 

Atmospheric radio noise is often tlie limit in^^ factor 
in radio coiuniunicatioiis at I'requeucies up to about 
30 Mc/s. In the design of coinnuiriications circuits, 
it is necessary to know the detailed statistical 
characteristics of interfering noise in order to cal- 
culate the clianuel capacity or error rate for a noisy 
channel. The cunuilative amplitude-probability 
distribution (APD) has been found to be a useful 
tool in such analyses [Montgomery, 1954; Watt, 
et al., 1958; Crichlow and Disney, to be published]. 

Atmospheric radio noise is a nonstationary ran- 
dom process whose cliaracteristics change not oidy 
with time but also with bandwidth. The APD is 
usually measiu'ed as a time sequence of several 
simultaneous levels, and the necessarily long integra- 
tion times make difficult the obtaining of a contin- 
uous curve. The need to overcome tliis difficulty 
has led to the developnu^nt of a method of predicting 
the APD from tlie three statistical moments of 
atmosplieric radio noise measured by the ARN-2 
noise recorder of the National Bureau of Standards 
[Crichlow et al., 1960a; Crichlow et al., 1960b ; Fulton, 
1961]. These moments are measured on a world- 
wide basis [Crichlow, 1957], and the data from these 
measurements are summarized and published quar- 
terly [Crichlow, Disney, and Jenkins, 1957-1961]. 

The three moments are measured for a power band- 
width of about 200 c/s, so the APD derived from 
them is valid only for this bandwidth. For this 
reason, a method for converting the APD of atmos- 
pheric radio noise at a 200 c/s bandwidth to a range 
of other bandwidths was developed. 

2. Parameter Definitions and Basic 
Assumptions 

It has been shown tliat the distribution function 
for atmospheric radio noise can be determined from 
its following three statistical moments [Crichlow 



et ah, 1960a]: 



rms voltage: r,.,!, 



average* voltage: iK 



1/2 



log of the volt agv : log Vu,^= .' I ^" log />/;>(/, t) 

where v is 1/V2 times tlie instantaiu'oiis envelope 
voltage, /; is the ])robability of v being exceeded, 
;/2— /i = 6 is the bandwidtli, and ^2 — ^i is tlie time 
interval of measurement. 

The three moments as measured are expressed, 
respectively as : 

/a=tbe effective noise figure 

= the external noise power available from an 
equivalent, short, lossless, vertical antenna 
in decibels above kth (the thermal noise 
power in a passive resistance at room 
temperature, t, in a bandwidth, 6; where k 
is Boltzmans' constant). 
V(i=Vs,we in decibels below Vrms- 



L, 



^log 



in decibles below Vr, 



Fa is independent of bandwidth for a uniform 
spectrum (a condition closely approximated in 
normal communications bandwidths), but Vrms varies 
as the square root of bandwidth. Va and La are 
always positive, since for this type of distribution 
function, ?^rm5> ?^ave>^iog. Becausc the shape of 
these distribution curves is dependent only on Va 
and L^*, which have been normalized to F^, it is 
possible to construct distribution curves of a giv^en 
form factor for various combinations of Va and La. 
Measurements of Va and L^a tend to sliow that 
La is a linear function of Va. Figure 1 shows La 
verstis Va for variotis bandwidths, seasons, and time 



* Unless otherwise stated, the moments used in the analysis are assumed to be 
the true moments of the distribution function. 
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Figure 1. L^ versus Vj from measured distributions. 



of day from measured distributions. The relation 
id=1.69 y^+0.72 seems to hold in gerenal. Because 
of the correlation between V^ and L^, one can deter- 
mine the most probable APD^s of atmospheric 
radio noise for the range of Va normally encountered. 
A family of these is presented as figures 2 and 3. 

The APD of atmospheric radio noise can be closely 
represented by a three-section curve on Rayleigh 
graph paper, i.e., two intersecting straight lines 
joined tangentially by a circular arc, figure 4. This 
particular Rayleigh graph paper has scales chosen 
so that a Rayleigh distribution function plots as a 
straight line of slope — K- These scales are labeled 
noise level in decibles above ^^rms versus the percentage 
of time that each level is exceeded. 

The lower section of the curve (low levels exceeded 
with high probabilities) represents the part of the 
noise composed of many random overlapping pulses, 
and plots as a straight hne Rayleigh distribution 
[Crichlow et al, 1960a]. The upper straight hne 



section of slope less than — }^ (high levels exceeded 
with low probabilities) represents the part of the 
noise composed of large, infrequent nonoverlapping 
pulses. If the exponent in the expression for the 
Rayleigh distribution function is raised to the power 
— jis, where s is the slope of this upper straight Hne, 
an expression of the upper portion of the APD will 
be obtained. For this reason, the upper portion is 
sometimes called a power Rayleigh. 

Four parameters are used to define the distri- 
bution as shown in figure 4. A is the decibel dif- 
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PERCENTAGE OF TIME LEVEL IS EXCEEDED 

Figure 2. Most probable amplitude-probability distribution 
functions for atmospheric radio noise. 

ference between the z^rms level and the Rayleigh line 
at 0.5 probability, and determines the amplitude of 
the Rayleigh section of the distribution. C and 
X fix the amplitude and slope of the power Rayleigh. 
C is the decibel difference between the Rayleigh line 
and the power Rayleigh line at 0.01 probability and 
X is the ratio of the slope of the power Rayleigh 
relative to the slope of the Rayleigh, that is, X= —2s, 
The parameter B describes the circular arc and is 
defined as the decibel difference between the inter- 
section of the two straight lines and the line tangent 
to the circular arc at its center. For atmospheric 
radio noise, an experimental correlation between B 
and X has been found, with ^=1.5 (X-1). These 
parameters have been determined as functions of 
Va and i^. 

The following assumptions are used to determine 
the relationship between the APD's for a sample of 
atmospheric radio noise received through different 
band widths. Each of these will be discussed in 
detail as they arise. 

1. The shape of the distribution for the prob- 
abilities of interest (10~^ and greater) will be of the 
above form for any bandwidth considered. 

2. The rms value of the distribution will vary as 
the square root of the power bandwidth, increasing 
with increasins: bandwidth. 
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Figure 3. Most probable amplitude-probability distribution 
functions for atmospheric radio noise. 

3. As the bandwidtli decreases and becomes quite 
small, the APD approaches a Rayleigh distri})ution. 

Since with the above assumptions the cunmlative 
distribution is determined by its moments {Va and 
Ld), the manner in wliich the moments vary with 
bandwidth needs to be determined. With Va^ and 
Zdj designating; the moments of the original APD at 
bandwidth 61, and Vu^ and L^^ designating the mo- 
ments of the desired distribution at bandwidtli 62- 

1. X2 and C2 nmst be determined as functions of 
Xi, Ci, and w, where w=^b2/bi and, 

2. Vd2 ^iid Ld^ must be determined as functions 
of Vd^, Ld^, and w from (1). 

3. Transformation of the Power Rayleigh 

Section of the Distribution 

The high-amphtude, low-probability section of the 
APD represents a train of nonoverlapping pulses. 
If the response of a receiver with bandwidth 61 to 
one of these pulses is a pulse of amplitude a and time 
duration t, and then the bandwidth is changed to 62, 
the response will be a pulse of amplitude wa and 
time duration t/w. For tliis section of the APD 
Fulton [1961] has shown that every point (p, v) 
corresponding to a bandwidth hx transforms to the 
point {pjw, wv) for a bandwidth 62, ^s long as the 
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Figure 4. Definition of parameters for amplitude-probability 

distribution for atmospheric radio noise. 
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receiver response maintains the criterion of nonover- 
lapping pulses. However, if the transformation 
{jp, v)^(p/w, wv) is made point by point, the results 
will not be another straight line, but a curve. As 
shown by the dashed lines in figure 5, the curvature 
is very slight in the probability range for which the 
power Rayleigh holds, because in the region in 
which the curvature becomes objectionable the 
criterion of nonoverlapping pulses is no longer met. 

Since it is desired to portray the transformed dis- 
tribution with the same form factor as the original, 
the curve will be approximated by an appropriate 
tangent to the curve. The slope and point of tan- 
gency of this line will be determined as a function 
of the bandwidth ratio w, so that a suitable approxi- 
mation is obtained; and a reciprocal relationship 
results, i.e., transforming from 61 to 1)2, and then from 
62 to 61, the original distribution is regained. 

The equation of a cumulative Rayleigh distribu- 
tion is: 



or 



p=exp 



1 



[-a 



logi>=-- [-log ( 



-lrii>)] + 2log^'^ 



where v is 1/V2 times the instantaneous envelope volt- 
age, J) is the probability of v being exceeded, and 
v^ is the mean square voltage. 

Since the slope of the power Rayleigh line is — —^ 

the equation of the distribution function for the 
power Rayleigh will be 



l9g^;=-— [-log (- 



1 , —. 



-In^?) 1+2 log 



2 



where Vt^ is the mean square voltage of the corre- 
sponding Rayleigh distribution, 



or 






Transforming the power Rayleigh by the band- 
width ratio w, 



or 



t(;^=exp 



X 






bg v=— log ( 



-\\\wj)) +- log wh?. 



Putting tlie above transformed curve in terms of 
our coordinates with 

x=— log (— Inp) 

or p=exp [—10"''], 

and ^=log V, 
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Figure 6. Ratio of X2 to Xi versus bandwidth factor w. 



it is seen that 



y=j\og (10- 



dy 
dx 



-u 



-Inw) +0 log w\^- 



Tn^ 



hip^lnw 



If Si{si = — Xi/2) is the slope of the original power 
Rayleigli and 82(82= dy/dx) is the slope of the result- 
ing power Rayleigh, 

X2_ \np 
Xi Inpw 

The point for the best tangent approximation will 
change with each transformation w. The point of 
tangency must be chosen such that the above men- 
tioned reciprocal relationship is maintained, and 
such that the transformed curve is well approxi- 
mated by its tangent in the interval of interest. 
This tangent point for each transformation was de- 
termined to be p^O.OOOl for wy>l and 61 = 200 c/s. 
The above requirements will be met if the tangent 
point is chosen at ^ = 0.0001 for 'U;>1, and the ap- 
propriate value for w<^l to maintain the reciprocal 
relation. 

In order to transform the power Rayleigh, then, 
the point (0.01^^;%, v) is translated to the point 
(0.01%, wv) and through this point the desired power 
Rayleigh with a slope given by figure 6 is drawn 
(w'^l). Figure 6 shows the ratio of X2 and Xi as a 
function of w. 

4. Transformation of the Remaining Section 
of the Distribution 

Having obtained X2 as a function of Xi and w, 
C2 must be determined as a function of Ci, Xi, and 
w. This is obtained by choosing various Xi and Ci 
for a particular w, transforming the power Rayleigh 
portion according to the above, and then determining 
O2 (that is, locate the Rayleigh portion of the desired 
distribution) such that the i^ms values of the original 
distribution and the desired distribution are related 
b^^ the V'^- An example of the procedure is given 
in figures 7 and 8. In figure 7, the distribution to be 
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transformed is chosen {Xi = Q, Ci = lO db), and from 
the known relationsliip between A, C, and X (fig. 9), 
Ai is seen to be 9.9 db and the relative ^rms level is 
arbitrarily set at the 20 db level. 

Transforming by the bandwidth factor 20, the 
point (Pi, Vi) is transformed to the point (Po, V2), 
where P2 is 0.01 percent, Pi is 20 (0.01%) or 0.2 per- 
cent, Fi is 36.5 db, and V2 is Fi + 20 log 20 or 62.5 
db. X2 is found to be 8.88 from figm-e 6 and the 
power Rayleigh section is transformed. The Frms 
level of the desired distribution transforms as the 
V20 or 13 db, so Frms2 ^^'^ ^^ '^*^ ^^- Using figure 9 
and choosing various C2, the correct 62 is found to 
be 9.6 db (fig. 8). Figure 10 shows C2 as a function 
of Ci and Xi for a bandwichh factor of 0.5. 

From previous work [Criclilow et al., 1960a] and 
the above, the following functional relations are now 
known: 

X=MVa,La) 
C =MV„ L,) 
X,=MXu w) 

These may be manipulated graphically to obtain 
Va,^ as a function of V^^ and w. The results are shown 
in figures 11 and 12. 

In order to deterinine the shape and amplitude of 
the APD for almospluM'ic radio noise at bandwidth 
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62 from Frmsi ^^^ ^^1 i^^easured at bandwidth 61, then, 
Fd2 is obtained from Va^ and w from figures 11 or 12, 
and the shape of the APD for this F^^ is found from 
figures 2 and 3. The amplitude^ of theAPD at band- 
width 62 is determined from the relation Frmso== 

Frmsi+lO log W. 
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5. Accuracy Considerations 

As the bandwidth is increased indefinitely, the 
above analysis becomes inaccurate since the tangent 
approximation is no longer sufficient. However, since 
Vd increases with bandwidth, and for large V^ the 
integrals defining Va are determined almost com- 
pletely by the power Rayleigh portion of the cumula- 
tive distribution ; the behavior of Va for large band- 
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Figure 10. C2 as a function of Ci and Xi for w=O.S. 
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18 



width ratios may be determined by investigating the 
behavior of the power Rayleigh. 

Since the power Rayleigii is coi^iposed of nonover-j 
lapping pulses, the character of the repsonse of ai 
bandwidth -limited circuit to this train of pulses can! 
be determined from the response to an individual! 
pulse. Therefore, the average value is independent! 
of bandwidth, and the rms value varies as the square) 
root of the bandwidth ratio. V4 (the deviation ini 
db between the average and rms values) will vary asi 
the square root of the bandwidth ratio, or Va will! 
increase 10 db per decade of bandwidth ratio w. The; 
above may also be shown by actual integration of the 
power Rayleigh and the transformed power Rayleigh. 

The above also gives a good check on the accuracy 
of the previous analysis. The results must approach 
this 10 db per decade law as w increases and also asj 
F^j increases. Figure 12 shows that this is indeed | 
the case. 

For a bandwidth of 200 c/s, it has been found! 
experimentally that the distribution saturates at a 
probability of 10"^, i.e., there will be almost no pulses 
with amplitude greater than that for which the 
probability of being exceeded is 10"^ This is not 
true for other bandwidths. In general, the distribu- 
tion will saturate at slightly lower probabilities for 
larger bandwidths and slightly higher probabilities 
for smaller bandwidths. By assuming the same form ^ 
of the distribution for all bandwidths (saturation at! 
p = 10~^)j an error is introduced. This error isi 
indicated by the shaded portions of figure 13. 

An indication of the maximum error in V^ can be 
obtained by evaluating the average and rms values 
for a power Rayleigh in the probability ranges [0,1] 
(that is, assuming no saturation) and [10~^,1]. This 
will determine the maximum possible error in the 
average value and the rms value. Since these two 
errors will be in the same direction, the diflference 
between them (in db) will be the error in Va. 
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Figure 12. Effect of bandwidth on Vj. 
where 
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Figure 13. Saturation of amplitude-probahility distributions 
for various bandwidths. 

We must evaluate 

^ave= J^^ vdp and V rL= J^^ vHp 



^=exp 



\r(0] 



and pi assumes tlic two values and 10 ^ Under the 
change of variable 



For :pi = 0, then 

For ^1= 10"^ ^'ave cannot be easily evaluated in closed 
form in terms of X. If X is chosen corresponding to 
the power Rayleigh with the largest slope likely to be 
encountered, ^ave can be evaluated. Choosing X=12 
or a slope of —6, 



i^ave 
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For the probability interval fO,l], and A''=12, 

Therefore, the maximum relative error in ^ave is 
0.13 db. 

In similar fashion, for the probability range [0, 1], 

?'rn.s=2.19X10^ V^- 
For the probabihty range [10"^ 1], and A^=12 

t— 12 19! -|13.82 

The maximum relative error in iv^is will, therefore, 
be about 2.07 db and the maximum possible error 
in Vd somewhat less than 2 db. The maxinmm 
error in La will be about the same, since the error 
in z^iog will be even less than that in v^y^. 

In practice, the error introduced will be much 
smaller than the 2 db, as it was assumed that the 
distribution never saturates in calculating the 2 db. 
Even so, an error of 2 db in Va for large Va will 
not noticeably change the shape of the cumulative 
distribution, but only the relative amplitude of the 
distribution. 

Since the amplitude of the distribution (zVms) is 
measured experimentally, we are only interested in 
the effect of the error on the shape of the distribu- 
tion; therefore, the assumption of a uniform point 
of saturation for all bandwidths is a valid one, since 
the shape of the distribution for the probabihties of 
interest will not be changed. 

In using the above, it should be remembered that 
the transformation depends on La being a hnear 
function of K. Figure 1 shows this is generally 
correct, however, there is some variation of La for 
a given Va, and this variation increases as Va in- 
creases. That is, in any particular case, it is pos- 
sible for the true La^ and the assumed true La^ 
(from fig. 1) to be significantly different, especially 
for large Va^. Also, Va^ itself is subject to error in 
measurement, and this error will be propagated by 
the transformation. The possible error in Va^ for 
a given Va^ subject to error can be quickly deter- 
mined from the transformation curves. For ex- 
ample, suppose Vai was measured to be 5 db subject 
to a possible error of ±K db (4.5<Ki<5.5) and 
it is desired to transform by w=\{)0. Va^ of 5 
transforms to Va.^ of 22.7, Va^ of 4.5 goes to Va^ 
of 21.7, and Va, of 5.5 to Va^ of 23.7. In this case, 
then, the i K db error in measurement has trans- 
formed to a possible ± 1 db error in Va^ (21.7 <iVa^ 
<23.7). On the other hand, if we transform the 
above Va, by w=Q.Q\, the ±K db error is seen to 
transform to a possible ±0.07 db error in V^^- ^^^' 



this reason, then, each transformation should be in- 
vestigated as to the effect of measurement error in 
Frij^ on Vd2 and the effect of this possible error in 
Va^ on the shape of the resulting APD. 

8. Conclusions 

Since measurements of the APD for atmospheric 
radio noise have show^n the assumed form factor to 
hold over a wide range of bandwidths [Watt and 
Alaxwell, 1957; Watt et al., 1958; Crichlow et al., 
1960a], the above transformation will usually give 
good results. However, in using the above, it sliould 
be remembered that an error in measurement of 
Frfj^ will generally propagate as a larger error in 
Va.2 as the bandwidth is increased and a smaller 
error in Va2 as the bandwidth is decreased. Also, 
the linear relation between Va and La is only a good 
approximation and subject to possible error, espe- 
cially for large Va. If a F^^ and La, are encountered 
significantly different from the pairs corresponding 
to the APD's of figures 2 and 3, the APD for this 
pair of Va, and La, may be obtained from NBS 
Monograph 23 [Crichlow et al, 1960b]. This APD 
may then be transformed by the method outlined 
above. 

Experimental verification of this method of pre- 
dicting the changes in Va with bandwidth has been 
obtained, but with some doubts as to accuracy since 
the noise changed character in the period covering 
successive bandwidth measurements. 

An improved method of measurement is currently 
under development using simultaneous magnetic tape 
recordings of the noise through different bandwidths. 
Results of the new measurements will be published 
upon completion. 
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